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Abstract 

We give evidence, by use of the Thermodynamic Bethe Ansatz approach, 
of the existence of both massive and massless behaviours for the (f)2,i per- 
turbation of the M3 5 non-unitary minimal model, thus resolving apparent 
contradictions in the previous literature. The two behaviours correspond to 
changing the perturbing bare coupling constant from real values to imagi- 
nary ones. Generalizations of this picture to the whole class of non- unitary 
minimal models Mp^2p±i) perturbed by their least relevant operator lead to a 
cascade of flows similar to that of unitary minimal models perturbed by 0i,3. 
Various aspects and generalizations of this phenomenon and the links with 
the Izergin-Korepin model are discussed. 
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1 Introduction 

The astonishing progresses of the last decade in two-dimensional confornial Quan- 
tum Field Theory (CFT), initiated with the work of Belavin, Polyakov and Zamolod- 
ckikov ||I[], have developed in many directions. One of the most interesting ones is 
the investigation of integrable perturbations of CFT f^, where a lot of amusing 
non-perturbative phenomena can be observed and the renormalization group flows 
can be reconstructed exactly in many cases from the S-matrix of the corresponding 
Factorized Scattering Theory. 

The bare action ^ of a CFT Ai perturbed by one of its relevant operators $ is 
given by 

A = AM + >^fd'^x^{x) (1) 

where Am is the (formal) action of the CFT. In what follows we shall use the 
shorthand notation TW + $ to refer to such a model. The perturbing parameter A, 
together with the parity of the perturbing operator, plays a central role in deter- 
mining the possible behaviours (massive or massless) of the integrable theory. For 
example, in the celebrated (pi^s perturbations of unitary minimal models (Mp-|-0i^3), 
the perturbing operator is parity even and the model is not necessarily equal for 
positive or negative A. Indeed, it shows a massive behaviour for negative A and a 
massless one for positive A, which have been widely studied in the literature. There 
are other models, like the unitary and non-unitary minimal models Mp^q perturbed 
by their 0i^2 or 02, i operators, as well as many other theories, where the perturbing 
operator is odd, and this implies an invariance for A —* —A telling us that whatever 
the behaviour for positive A is, it will be reproduced for negative A too. 

The interest of this paper is to clarify when a certain integrable perturbation of a 
CFT admits both massive and massless behaviours. The importance of discovering 
new massless behaviours lies in the fact that they interpolate between non-trivial 
ultraviolet (UV) and infrared (IR) CFT's, thus providing very interesting infor- 
mation on the structure of the two-dimensional Renormalization Group Space of 



Actions. It has been widely believed that the two behaviours can coexist only 
when there is no symmetry A -^ —A, i.e. when the perturbing operator is even. 
When this symmetry is present, instead, one could conclude at first glance that 
only one of the two possible behaviours is allowed. However, Fendley, Saleur and 
Al.Zamolodchikov [P, Q have recently given support to the possibility of a massless 
flow in the Sine-Gordon model where the perturbing operator is Z2-odd. This new 
type of flow is related to an imaginary coupling constant in the potential. In this 
paper, by examining the simple M3 5 + 02, 1 theory, we shall learn that such a subtle 
situation can occur also in minimal models and other rational CFT's, thus enlarging 
considerably the set of possible integrable models. As a by product, this will also 
suggest the existence of massless flows with imaginary coupling constant (similar 
to those described in [^, Q for Sine-Gordon) in the Izergin-Korepin model. 

The plan of this paper is as follows: in sect. 2 we summarize some results about 
the Ms^s + 02,1 theory. Arguments in support of a massive behaviour for this theory 
have been given in ref. |J5|, while a massless interpretation has been supported in 0. 
The oddness of the 02, 1 operator seems to imply a contradiction between the results 
of [§ and [§. 

In sect. 3, starting with the hypothesis that the model is massive, with the S- 
matrix proposed in p, we give a set of Thermodynamic Bethe Ansatz (TBA) 
equations governing the Renormalization Group (RG) evolution of the Casimir en- 
ergy of the vacuum and of the first excited state on a cylinder. This is obtained by 
folding the known set of TBA equations of a suitable lys-minimal model perturbed 
by its (f)id,adj operator. The TBA set passes various checks that give a high level of 
confidence to our conjecture. 

Then, in sect. 4 we turn our attention to the other possibility, i.e. to the mass- 
less behaviour. We show that it is possible to construct a self-consistent massless 
scattering theory interpolating the ultraviolet (UV) M3 5 model and the infrared 
(IR) M2,5, the celebrated Lee- Yang singularity. As this scattering theory turns out 
to be diagonal, the TBA equations can be deduced easily in this case. They corre- 



spond exactly to the "inassless version" of the previous massive TBA, according to 
the empiric rule introduced by Zamolodchikov to pass from massive to massless 
TBA by substituting mass terms with left and right movers in a suitably symmetric 
way in the TBA equations. The perturbative analysis of the TBA solutions for the 
ground state Casimir energy leads to the solution of the original puzzle: in this case 
it is the X ^ iX transformation that leads from massive to massless regime. 

Sect. 5 deals with the generalization of the Ms^s + 02, i result to a whole series of 
models Mp 2p±i perturbed by their least relevant operator and the link of this series 

(2) 

of models with the Izergin-Korepin A2 affine Toda Field Theory. A straightforward 
generalization of the M3 5 TBA equations leads to the picture that for each of these 
models there exists both a massive behaviour, for which the S-matrix should be 
deducible as a folding of a suitable PFs-minimal model one, and a massless regime 
of interpolating flows, accompanied, like in the usual Mp + 0i^3 case, by a staircase 
that has already been reported in |P. We also briefly discuss why in this case, 
unlike the Mp + 0i 3 one, large p perturbation theory cannot be used to study the 
massless flow. 

Our conclusions on these results are collected in the final sect. 6, where we collect 
comments about other theories that could share the same phenomenon and try to 
propose some criterion to select possible candidates. Here we also briefly discuss 
some open issues and possible developments for further research. 

2 The M3 5 + (/)2,i puzzle 

We begin our investigation by presenting an intriguing puzzle that appears in the 
02,1 perturbation of the M3 5 non-unitary minimal model, one of the simplest CFT's 
after the Lee- Yang singularity (M2,5) and the critical Ising model (M34). The 
^3,5 CFT has central charge c = — | and three non-trivial scalar primary fields 
01,2 (a^), 01,3(3^) and 02,1 (a;) of (left) conformal dimensions Ai^2 = ~^) ^1,3 = | a^id 
^2,1 = |, which satisfy the well known fusion rules of minimal models. The ground 



state of the theory does not coincide with the conformal vacuum |0) as in unitary 
models, but better with the state \Q) = 0i^2(O)|O). The ground state Casimir energy 
on a cyhnder is therefore not proportional to the central charge c, but to the so 
called effective central charge c = c — 24Ai 2 = |. 

In the following we shall consider the model M^^^ + (/)2,i, with action 

A = A/3,5 + A /" rf2x02,i (x) (2) 

that can be interpreted in Statistical Mechanics as the thermal perturbation of the 
Q = 4 cos^ ^ critical Potts model. A counting argument a la Zamolodchikov and 
an explicit construction of conserved currents for the simplest cases shows that this 
model is integrable. Moreover, the operator 02,i is odd, as one can see from the 
symmetries of the fusion rules and therefore the model is invariant for A -^ —A. 

A factorizable S-matrix for the model (0) can be proposed by quantum reduc- 
tion of the Izergin-Korepin model, as explained for general Mp^g + 02,i in ref. |§| 
and specialized to this simple case in ref. |^. The model presents two vacua (call 
them and 1). The asymptotic states are identified with kink states \Kqi) and 
|i^io) interpolating them and a single state over the 1 vacuum \Kii). There is no 
similar state |-ft"oo) over the vacuum instead. The admissibility diagram of the 
corresponding S-matrix is depicted in Fig. |1], while the detailed expressions of the 
S-matrix elements are given in ref. [^. The author of ref. ||^, after supporting his 
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Figure 1: The T2 diagram describes the vacuum structure of the M3 5 + 021 modeL 
hypothesis that the model is massive by a truncated conformal space analysis, con- 



eludes that this is the only possible behaviour of the model, due to the symmetry 

However, Ms^s + 02,i belongs to a peculiar series of models, that we shall denote 
below as Mp + ip, (p > 4) where 

, Mp+i „ for p odd 
Mp={ —" "^ (3) 

Mv^pj^i for p even 



and 



;/)2,i for p odd 
-0 = ■^ ' = least relevant operator of Mp (4) 

/)i^5 for p even 



The effective central charge of this series of models is given by 
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(5) 



p{p + 1) 

From the Statistical Mechanics point of view, this series describes a possible mul- 
ticritical generalization of the Lee- Yang model, definitely different from the one 
defined by a chain of usual 0i 3 perturbations and ^3 1 attractions.]^ 

In ref . |0 , evidence is given for a staircase model related to the A2 Affine Toda 
Field Theory, which flows close to all models Mp. It is then natural, by analogy 
with all the other known examples of staircase models [^ |T^, |rT| to conjecture that 
its 9q ^ 00 limit mimics a series of massless flows 



Mp + ■?/'—>• Mp_i attracted by < 



02,1 for p — 1 > 4 even 

01,5 for p — 1 odd (6) 

TT for p — 1 = 4 



Notice that the attracting operators are the least irrelevant operators in the Mp_i 
models. The situation looks very similar to the one known in unitary minimal 
models where the perturbation by the least relevant operator 0i,3 at ultraviolet 



^The interpretation of these niulticritical versions of the Lee- Yang model as some Lee- Yang like 
singularities related to each of the multicritical Ising points corresponding to the unitary minimal 
series Mp is not straightforward and needs further investigation. 



(UV) leads, for positive A, to an IR point attracting the flux by its least irrelevant 
operator 03^i. The considerations of ref . |l^ should apply here similarly. 

In particular, this supports the hypothesis of the existence of a massless be- 
haviour for the M3 5 + 02,1 model, flowing towards an IR limit governed by the 
M2,5 model (the Lee- Yang singularity) with attraction operator given by TT. (T 
and T denoting the two components of the stress-energy tensor of M2,5). The most 
dramatic point about this is that the invariance A -^ —A should then imply that 
the model is always massless, in striking contradiction with the results of ref. [Q . 

In conclusion, the M3 5 + 02, 1 puzzle can be summarized as follows: the A -^ —A 
invariance seems to predict only one possible behaviour (massive or massless) for 
the M3_5 + 02,i model. However, in the literature some reasonable evidence has been 
given to support both the massive or massless interpretation. Who is right? 
Or, is there a way out allowing both behaviours without contradiction? This is the 
main question we shall clarify in the present paper. 

3 TBA for massive Mgs + 02,i 

As we said in the previous section, a massive scattering theory has been proposed 
for the model in question in ref. |^. The corresponding factorized S-matrix is 
not diagonal and a direct Bethe Ansatz approach to deduce the TBA equations 
governing the evolutions of energy levels on a cylinder along the RG flow is very 
hard. One can however avoid this difficulty by guessing the TBA equations and 
then verifying if they pass various checks to give them credibility. The starting 
point of our guesswork is the observation that the central charge of the models 
Mp, eq.(^, is exactly half of that of the {W^)p series of minimal models of W^- 
invariant CFT [l^. This suggests that a relation could occur between the (^^3)^ 
models perturbed by their (t)id,adj operator, {Wz)p + (l)id,adj for short, and our Mp + ip 
integrable models. The case of the bottom models of the two series, the scaling 
Potts model ((14^3)4 + 0a=2/5) and the scaling Lee- Yang singularity (M2,5 + ^1,2) 
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has been dealt with in ref. [0, where it is stressed that the folding procedure 
leading from the Potts S-matrix to the Lee- Yang one induces an analogous folding 
in the structure of TBA equations. The second models in the two series are the 
(^^3)5 + 4>id,adj and the M3 5 + 02, 1 models respectively. The S-matrix of the former 
is given in [|l^ and one can verify against the explicit expressions of ref. P] that 
the same kind of folding procedure works here too to produce the S-matrix of the 
latter. It seems then natural to conceive that the TBA equations for the Ms^s -|- 02,i 
model are those of the (1^3)5 + (pid,adj suitably folded. 

These latter are given in [jl6|, |l^] and read as follows (here -k means convolution 
and A^ = log(l + e^l), L{ = log(l + e"^'-)) 



Zn 



EG 

6=1 



ab/ 



[K 






1,2 



1,2 (7) 



where the kernel d) has the form 



2 cosh p 



the two matrices G and H have non-negative integer entries (and therefore can be 
thought as incidence matrices of two graphs Q and H respectively) 



G = H 



1 

1 



n = A2 



(9) 



and the energy terms I'aiO) are chosen as follows 
ul{e)=mR cosh e ul = 



1,2 



(10) 



where m is the mass of the fundamental kink of rapidity 6 and R the radius of the 
cylinder on which the theory is put. 

This complicated set of equations can be more clearly depicted in graphical 
form, as usually done by many authors (for details see e.g ref. [l^) by drawing the 
tensor product graph Q x H (A2 x A2 in our case), whose connectivity reproduces 
the coupling of variables in the TBA eq.(|^ and giving to each node {i, a) a choice 
for the energy term z/*. Our case is depicted in Fig. ^ja. 




Z2-folding 



(b) 



Figure 2: Graphical representation of TBA systems encoded on Q xli.. The diagram Q develops 
vertically, H horizontally: (a) the A2 x A2 case, corresponding to (^3)5 + (f>id,adj', (b) the Ti x A2 
case, corresponding to Afa^s + (f>2,i. Nodes denoted by O are attached an energy term 1^1 = 
(magnonic nodes), those denoted by • are attached i^l — mR cosh 9 (particle nodes). 



The TBA eq.(^ is symmetric in the a indices. The aforementioned folding 
procedure gives our conjecture for the TBA corresponding to the S-matrix of ref. |^ 



£* + — 0* 

Zn 






1,2 



:ii) 



where u^ = mRcosh.6 and u'^ = 0. It can be interpreted in graphical form as the 
folding of the Q = A2 diagram by its Z2 symmetry to the Ti diagram. Then this 
TBA set can be depicted on the Ti x A2 graph (see Fig. ^). 

Now that we have our TBA equations at hand, we can evaluate various quantities 
with them to learn more about the theory described by the S-matrix of [^]. The 
most immediate thing to compute is the scaling function 

3 



c[r] 



- / deu\e)L\e) 

TT^ J-00 



(12) 



where r = mR is a dimensionless scale, such that t = log r can be interpreted as the 
RG time going from —00 at UV to +00 at IR. In particular the UV limit r ^ can 
be exactly computed by resorting to Dilogarithm identities, and results in cuv = § 
as expected. Moreover, each solution of the TBA equations ( pA]) is also a solution 
of the following system of functional equations 



Y' e 



m 



Y'ie 



y 



i + Y\e) 



n,ti(i + r^W-^) 



l\W:> 



(13) 



The following periodicity property is shown by this system 

Y\0 + 2'Ki) = Y\e) (14) 

Al. Zamolodchikov [|1^ was able to relate this periodicity to the conformal dimen- 
sion of the UV perturbing operator. In our case this turns out to give exactly 
^pert = 4 as expected. These two fundamental signals {cuv = f and Apert = |) 
give us confidence in our conjecture, but many other interesting properties can be 
analyzed. We try to summarize them in the following. 

3.1 Comparison with perturbation theory 

Around UV, the behaviour of the scaling function c{r) is well known to be summa- 
rized by the formula 

oo 

c{r) = CUV + bulk + ^ Cn{rT (15) 

n=l 

where y = 2 — 2Aperi = |- The non-perturbative bulk term that takes into ac- 
count long range fluctuations can be calculated from general principles in the TBA 
approach []T1|, ^, ^, |^, |2^ and in the present case takes the form 



3 

bulk = — -r^logr (16) 

47r^ 

The scaling function c{r) can be computed numerically up to very high precision. 
Once the bulk term is subtracted, the c„ coefficients can be estimated via polyno- 
mial fit and compared against the results of perturbation theory around the UV 
conformal point. The UV perturbative series is 

oo 

dpert = CuV+T.PniXR'^^r (17) 

n=l 

with the coefficients P„ given by the UV correlation functions 

( — l)"" i-n-l j2^ n-l 

^n = 12^v^y n ^2vr|z.|)i/2 <^l'^^-^^^'^)n'/'2,i(^„%)|^) (IS) 



A (2^l^^-l)'^ 



Therefore the coefficients c„ can be compared against the Pn once the proportion- 
ahty constant k = Xjvn}!'^ is known 

Cn = t^'^Pn (19) 

Recently Fateev [^ used external magnetic field techniques to calculate k indepen- 
dently of TBA for a large class of theories. His result in our case can be summarized 

as 

r K i2 r^i2 

rr) i TTT J- ~ 

(20) 



K 



2 ^ 



r2. 



2 

r 



thus providing the proportionality constant between A and m}''^ 

A = 0.253001im^/^ (21) 

Surprisingly, for the first non-trivial coefficient P2 (the P„ with n odd are all 
zero by the symmetries of the fusion rules) the check can be done even if k was not 
known. This is due to the curious result that P2 = although the 4-point function 
needed to calculate it is not zero: 

/dp Z 
^^(01,2(00, OO)02,1(1, l)02,l(^, ^)01,2(O, 0)) (22) 

This 4-point function involves only one conformal block and its form can be fixed by 
using informations on its monodromy and the SL{2, C) invariance. As a final result 
one is led to compute a Dotsenko- Fateev like integral whose quadrature is identically 
zero thanks to the property of the hypergeometric function F{—l,b; —b; —1) = 0, 

What we have done was first to fit {c{r) — bulk) with a polynomial having both 
even and odd powers of r^/^, to check that really the c„ with n odd are all zero 
within numerical error. Then we turned to a fit with only even powers of r^/^ 
to make more precise estimates of the C2„. The results are collected in the first 
column of table |I], where it appears clearly that the C2 coefficient is zero within an 
approximation of 10"'^'^. Table ^ contains also other information that will become 
clear in the next pages. 

10 



Cn 


massive 


massless 


Ac 


Co 


0.6 


0.6 


/ 


C2 


-5.4757325311081032 • 10'^^ 


-8.8531745297863157-10-15 


±6 - 10-1^ 


C4 


-0.1755364844390 


-0.1755364844395 


±7-10-13 


c& 


2.2294559673 ■ 10"^ 


-2.22945596576 - 10-^ 


±4 - 10-12 


Cg 


-1.012291-10-3 


-1.012290- 10-3 


±2 - 10-9 


ClO 


-8.6210-10"^ 


8.6212-10-^ 


±2 - 10-s 


C12 


-9.87-10-5 


-9.85 - 10-5 


±3 - 10-^ 



Table 1: 

Massive and massless UV coefficients. The last column reports the difference between 
the absolute values of the two. 



3.2 TBA for the first excited state 

The massive TBA system in presence of a chemical potential (with suitable chioces 



of the latter) is known to yield the behaviour of various excited states pj, ^ . We 
find that the TBA system obtained from (|l^) with the substitution 



-e{0) 



-eie) 



(23) 



in A* (6*) and L\9), describes the first excited state of the perturbed massive the- 
ory. The UV limit of this state is the conformal vacuum |0), therefore its Casimir 
energy is directly proportional to the central charge of the M3 5 model. The sub- 
stitution ( P5D does not influence the periodicity property of e-^*-^-* and in the UV 
regime we have a perturbative expansion like eq. (|l^). At r ^ we obtain the cen- 
tral charge of the theory via standard dilogarithmic sum-rules. We find cuv = ~| 
as it must be for the model M35. We also solve numerically this TBA system and 
in Fig. ^ we report, as functions of R, the first two energy levels obtained from 
the numerical solution of eqs. (|11|) for the ground state energy Eo{R) and with the 
substitution (^3|) for the first excited state energy Ei{R). We see that the two levels 



11 



exponentially degenerate in agreement with the double well potential of the theory 
as interpreted in p. 

In Fig. ^ we represent the function AE{R) = Ei{R) — Eo{R) compared with the 
result obtained in ref. using the Truncated Conformal Space Approach (TCSA). 
Up to a little discrepancy due to the finite level truncation of the Hilbert space in 
TCSA, we see that the two results agree, thus giving another piece of strong support 
to the TBA system ([Tl|). 

In order to give even more support to the massive TBA we check the first nu- 
merical coefficient against perturbation theory along the same lines of the previous 
sub-section. In this case the value for k, is needed and gives 



sr 



2 



r 



2 



r 



ci = kPi = ^-^ ^-^ — ^ ~ 0.551329 (24) 



"l' 

.4. 


2 

r 


" 7 ■ 
.12. 


2 

r 


"5" 
.3. 



Using standard numerical solution of the excited TBA system we find ci = 0.5513 ± 
0.0001, in very good agreement with the perturbative result. 

3.3 IR regime 

For massive TBA systems, we may compare the IR behavior predicted by TBA 
equations and the one predicted by general considerations (cluster expansions) on 
the associated scattering theory. For kink like scattering, these latter predict that 
the leading contribution to E{R) is proportional to the integral 

/(r) = -— f coshee-'^^°''^^d9 (25) 

27r J 

multiplied by the largest eigenvalue A^ax of the kink adjacency matrix, as discussed 
in ref. [^ . So the leading term of the IR asymptotic is 

Eo{R) ~ Ama.ml{r) (26) 

The vacuum structure of our theory is described by the T2 diagram of Fig. |I|. Its 
incidence matrix has eigenvalues 

Ajnax = '2 cos (-] , A^j^ = 2 COS f — j (27) 

12 



E(R) 




0.5 1 1.5 2 2.D 3 



R 



Figure 3: The vacuum and the first excited states obtained using the TBA system. 
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AE{R) 




R 



Figure 4: The energy gap between the vacuum and the first excited state decaying exponentially 
with the volume R. We compare the gap shape obtained using the TBA with that obtained using 
TCSA (dotted fine). 
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Our TBA system (|Tl]) in this limit predicts 

e{6) r^ r cosh 9 — log [l + e~^^°°n = r cosh 9 — log 
and 



2 cos I — 



Eq{R) ^ f d9 cosh 9e-'^'^^ = 2cos (^) ml {7 



f28) 



(29) 



in agreement with the expected result. The lowest excited state should become 
degenerate with the ground state in infinite volume, from our TBA system modified 
by (p3D for the excited state we find 



e{9) ~ r cosh 9 — log (e~ 



■e(oo) 



r cosh 9 — log 



—2 cos 1 — 



(30) 



and 



-6(6) 



/'37r\ ^, , 
2 cos ( — ) ml [rj 



(31) 



777 r 

Ei{R) d9 cosh 9e- 

2tt J 

This result confirms that the splitting AE{R) decays exponentially with the volume 
R. Notice that the coefficient in ( pT]) corresponds to the second eigenvalue of the 
incidence matrix of T2. This seems to be a quite general property of the degenerate 
vacuum state theories and suggests the existence of some kind of general relation 



between TBA and its graph encoding [|T^ and the usual graph theory of S-matrix 
adjacency. This investigation is in progress. 



4 Massless Scattering for M3 5 + 



2,1 



M 



2,5 



In this section, we wish to show that it is possible to propose a massless scatter- 
ing theory describing a fiow from the UV M3 5 + 02,1 down to the IR M2,5 model. 
The S-matrix formalism can be developed for a massless scattering theory following 
the lines proposed by Alexander and Alexey Zamolodchikov in p6[ and developed 
in 1^, |2^. We refer the reader to those papers for a treatment of the subtelties 
concerning the definition of a massless S-matrix and the properties related to it. 
In particular the deduction of massless S-matrix for the minimal models described 
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in IQ will be followed quite closely. A massless scattering is defined as consis- 
tency requirement of the Bethe equations, thus overcoming kinematic difficulties. 
Instead of massive asymptotic states (particles), one introduces right (+) and left 
(— ) movers whose energy can be parametrized as ye^^ respectively (m is a scale 

of the theory). Four scattering matrices have to be considered: 5*++, S , 5+- and 

S-^. Parity invariance requires that 5++ = S and by analogy to the Mp + 0i^3 

case we also require 5'_+ = S'+„. 5++ is scale invariant and can be thought as the 
scattering matrix for the left sector of the IR CFT (M2,5 in our case). Analogy with 
the unitary minimal model series suggests to take it formally equal to the S-matrix 
of the IR model perturbed by its least relevant operator in the massive direction. 
Therefore, let us consider the massless scattering defined by taking 5++ = Sly, 
where Sly denotes the S-matrix of the massive perturbation of the M2,5 model by 
its operator of conformal dimension — | given in [^. i.e. let us choose 

We use here the notation (x) = '^■\ i fn^^-^^\ ■ ^J analogy with what assumed in 
ref. 0], we take S^^{6) proportional to 5++ (6' — ia), for some rapidity shift a, 
and S^^ oc S^^{6 + ia). A consistent choice, compatible with the requirement 
S-j^ = S^-, is to fix a = it/ 2 and the proportionality factor equal to 1. Therefore 

1\ / 2\ _ ,_i 



S^-{0) = - (-3 j (-3 j = (Sly)-' (33) 

We propose as scattering theory corresponding to the massless M^^^ + 02,i model 
the one governed by the 5*++ and 5+^ described above. 

4.1 TBA and UV behaviour 

Of course this proposal must be checked by computing the UV and IR behaviour of 
the theory defined by this S-matrix and showing that UV really gives Ma^s perturbed 
by 02,1 and IR gives M2,5 with attraction operator TT. This can be done by resorting 
again to the TBA approach. The deduction of the TBA system from the S-matrix 
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is particularly simple here, as the latter is diagonal. Standard calculations give the 
TBA equations for the scattering theory defined by (^|33D in the form 



^e±^ = e^{9) + ^[^^{L^ - L^)W) (34) 

where R is the radius of the cylinder and L± is short for log(l + e~^*). The 
convolution kernel 0(6') is given by 

m = -^-^\ogS^M = ^^\og 3^46) (35) 

Simple manipulations involving a Fourier transform, bring this TBA system into 
one having the same form of eg . (pUD , where now i = +,—. This massless TBA is 
therefore encodable, like the massive one of the previous section, on the Ti x A2 
graph. The only difference is in the choice of the energy terms attached to the two 
nodes. Here they read u^ = |exp(±^). The scaling function 

3 r 



c[r) 



±L fde[e'L4e) + e-'L4e)] (36) 

TT^ Z J 



can easily be evaluated in the two r — > (UV) and r ^ oo (IR) limits by resorting 
to Rogers Dilogarithm sum rules and gives cuv = | and cjr = |. The TBA system 
(0) can also be recast in the form of a set of functional equations equal to that of 
eq. (^) (we just pick up a different solution with another asymptotic condition on 
the pseudoenergies £„), thus predicting from the already mentioned periodicity that 
the UV perturbing operator has conformal dimension Apert = f- The stationary 
(i.e. independent of 6) version of the Y-system is a set of algebraic equations whose 
solutions, once inserted in Dilogarithm sum rules, not only predict the UV effective 
central charge, but also the conformal dimensions of other excited states and even 
subsets of the CFT fusion rules PD|. Such an exercise applied to the present case 



confirms with no doubt that the UV limit of the scattering theory ( p2| , p3|) is the 
M3 5 minimal CFT perturbed by its 02 1 relevant operator. 
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4.2 IR behaviour 

What is different from tlie previous massive case is tlie r ^ oo beliaviour near IR. 
Tfiis latter is dictated by tlie Y-system witli tlie right mover deleted, i.e. by the 
Ti X Ai system encoding the effective central charge, conformal dimensions and 



fusion rules |2^ of the c = | Lee- Yang singularity model M2,5. The Y-system peri- 
odicity at IR must be compared against the asymptotic expansion in j"-2{i+Aattr) j|^^ 
where Aattr is the conformal dimension of the irrelevant operator of the M2,5 CFT 
attracting the flow. It turns out that Aattr = 2, thus confirming the identification 
of such operator with TT, as expected. 

Results of perturbation theory around the IR Conformal point for the massless 
TBA (see ref. [^) are described by the asymptotic perturbative series 





Cpert ~ 


00 

n=0 




n 




(37) 


The first 6„ coefficients are 














60 = {l2Amin - Cir) 


bi = - 


{12Ajnin - Cmf 


, h-- 


= 2(12A™,„ - 


cm)' 














(38) 


We must compare the perturbative 


expansion (p^) 


with the numerical solut 


ion of 


our TBA system 
















Cpert > ; 

^ — n 


hn{RM)- 


n 




(39) 


The first three coefficients 


are 


&o = 


2 

5 






(40) 


(obtained using the dilogarithm sum-rule) 


and 









61 = 0.58041579 ± 8 ■ 10"*^ , 62 = 1-68440 ± 2 ■ 10"^ (41) 

Setting Kg = in~^ we use the second coefficient in order to fix the constant K ,we 
find K ='^hi = 3.627598 ± 2 ■ 10"*^ thus allowing to check the third coefficient. We 
find 62 = 2^^ [Ij = 1.68441 ±2 ■ 10"^ in good numerical agreement with (|^). All 
these checks confirm that the IR limit of the scattering theory (^,^) is the M2,5 
Lee- Yang CFT, with TT attraction operator. 
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4.3 Analytic continuation and solution of the puzzle 

Having now a quite solid ground of evidence for both massive and massless regimes, 
we turn our attention to the main question of the paper: how can the two behaviours 
coexist, in spite of the A —>■ —A invariance? To give an answer, first of all notice that 
the TEA systems ([Tl|) and ( ^Ij) are a perfect example of the recipe observed in a 
lot of integrable perturbed CFT having both massive and massless regimes, where 
one passes from massless to massive TEA by keeping the same structure of the 
equations, and modifying the energy terms only. Instead of the right mover |e^ one 
puts a massive particle energy term r cosh^, and instead of the left mover |e~^ one 
replaces a magnonic term with energy 0. This procedure is well known in unitary 
minimal models perturbed by (pi^s pT] , 0, as well as for many coset generalizations 
of them 122], |1^ and also on other examples ||3y, ^ |T^. This observation is also 



expected to be equivalent to relate the massive and massless scaling functions by 
analytic continuation in the parameter A [0]. For the Triciritical Ising model 



examined in [21], as well as for all Mp + 0i 3, the analytic continuation was from 

positive to negative A. Here we know this cannot be the case, due to the A -^ —A 

invariance. To have an idea of what kind of analytic continuation is needed, we just 

examine the perturbation theory around UV. For the massive case this has already 

been done in section 3.1, for the massless case the calculations repeat exactly the 

same structure. The non-perturbative bulk term, due to the different asymptotics 

chosen on the TEA diagram, is not the same as in the massive case, but rather 

reads 

bulk = — — + --r^ log r (42) 

47r Att^ ^ ^ 

Notice that in addition to the logarithmic contribution of the massive case we 
have here also a contribution proportional to r^, more similar to those discussed 
MM- 



m 



Once this bulk term is subtracted, the numerical criteria as in the massive case 
can be adopted here to extract the coefficients c^, (see eq.(0)) listed in the second 
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column of table |1[ Observe that, within numerical error, one passes from massive 
c{r) to massless one just by readjusting the coefficients as 

Cn ^ i-ir/'Cn (43) 

The Pn coefficients computed from CFT correlators must be strictly the same, which 
means, in view of the relation ([T9|) between c„ and P„, having the following recipe 
in passing from massive to massless regime 

X^iX (44) 

The scaling function (and therefore the ground state Casimir energy) remain real 
after the substitution (0), as it is (at least perturbatively) an even function of A, i.e. 
depends only on A^. Thus it is natural to propose ( PD as the analytic continuation 
allowing to pass from massive to massless regime. Notice that in view of eq. (pID , 
it is the massless theory that has real A, while the massive flow develops on the 
imaginary A axis. 

The analytic continuation now proposed can be checked even beyond the conver- 
gence radius of perturbation theory by resorting to the numerical technique of Fade 
approximants. Although we have only done some rough checks with this method, 
they agree, within numerical error, with the data for c{r) extracted from the two 
TBA's. 

The solution of the puzzle of sect. 2 is then that both papers [Q and [^] report 
part of the truth, in the sense that both behaviours are allowed from a physical 
point of view. For both of them it is possible to propose a scattering theory and a 
TBA system driving the RG behaviour from UV to IR. What is new in the point of 
view of this paper is to allow the perturbing parameter A to take in principle any 
complex value (which does not affect the integrability properties of the model) and 
then check if there is some direction in the complex A plane where one can still deflne 
a consistent theory having real energies for all the states. In the case M3 5 + 02, 1 
examined here we observe that the Z2 symmetry implemented by A ^ — A implies 
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that all observables are even or odd functions of A. Moreover, in a non-unitary 
theory like Ms^s + 02, i it is not necessary that the action (|^) must be real (which, as 
the operator 02, i is self-conjugate, would have implied reality of A). For example, the 
transformation X ^ iX keeps all even observables unchanged and trivially multiplies 
all odd observables by an i prefactor that can be readsorbed in the normalization. 
The theory defined by iX is then still a consistent one. 

We think that the behaviour observed in this simple theory is just an example 
of a phenomenon with a much wider scope, that can enlarge the zoo of integrable 
models by a lot of interesting new cases. This is somewhat parallel to the direction 
of investigation pointed out recently by ref . |^, ^ . In next sections we start exploring 
this larger zoo by discussing some other models sharing this phenomenon. 

5 The series Mp+ip and the Izergin-Korepin model 

In this section we generalize the results of Sect. 3 and 4 to the whole class of models 
M+ip. First of all, let us observe that all these models can be seen as quantum group 
reductions of the Izergin-Korepin model (the A2 Afiine Toda Field Theory) |§]. 
The latter has Lagrangian 

L = lid,<Pr + ^,e^''-^,e-^-' (45) 

and can be seen as a complex Liouville theory perturbed by the field V^cj)) = e~*2'"^. 
After quantum group reduction, we can consider the system ( ^5]) at £j: = ^ as a 
perturbed minimal theory Mp.g. In this picture the field e~^2't' is the perturbing 0i^2 
operator. At quantum level, by interchanging the role of p and q the field V{(j)) can 
also be associated to the 02, 1 operator. The theory (|i5| ) can describe both these two 
different integrable perturbations. This seems to provide only half of the models 
Mp + -0, those where ip = 02, 1. However also the other ones with ip = 0i 5 can be 
described in the same framework. One just has to change the role of the two vertex 
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operators in ( ^Sj) and send 13 ^ 213 to have 

L = \{d,<t^f-j,e-''' + j,e^''' (46) 

With this simple modification at ^ = ^ we can associate the Lagrangian (03) to 



Stt q 

the perturbation of the minimal theory Mp ^ by the operator 0i5. This operator, as 
the 021 operator, is relevant only in particular minimal theories, systematic check 
using the counting argument reveals that the field 0i5 is always integrable but it is 
relevant only for 2p < q, condition which is always satisfied by the models Mp with 
p even. 

We prefer to keep the interesting problem to write the (massive and massless) S- 
matrices for this series of models out of the scope of this paper. Here we just observe 
that they could be deduced both as reductions of the Izergin-Korepin S-matrix (the 
problems with real analyticity in ||^ seem to be circumvented in ref. [§^), or as 
foldings of the VFs-minimal model ones. The matching between the two formulations 
of the problem should give an interesting check on the ideas of this section. 

A simple proposal for a TEA system for these models can be made by just 
adding more and more magnonic nodes to the TEA of the Ma^s + 02,i model. This 
is compatible with the folding of W3 minimal models mentioned in sect. 3, and with 
the procedure usually adopted in minimal model series pT| , |^ as well as in many 
series of rational CFT perturbed by (f)id,adj operators [|l^ . The TEA system for the 
model Mp + ■?/; is written exactly as eq.(|ll|), with the sum now running from 1 to 
k = p — 3 and can be encoded on a Ti x Ak Dynkin diagram. The massive regime 
will be described by the chioce of energy terms with the first node as z/^ = r cosh 6 
and all the others i/* = 0. The massless behaviour, instead, is reproduced by 
choosing u^ = |e^, u'^ = |e~^ and all others z/* = 0. This is very similar to the 
TEA structure for unitary minimal models perturbed by 0i,3 introduced in |21, |^. 



Of course Dilogarithm sum rules predict the expected UV and IR effective central 
charges and the Y-system periodicities are compatible with the conformal dimension 
of the perturbing operators as well as the attracting ones. In particular the massless 
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regime consists in a series of flows hopping from Mp to Mp_i, exactly as described 
in Sect. 2, eq.(0). As in all known cases of series of hopping flows, it is not surprising 
that this is accompanied by a staircase model, which is the one described in ref. 0. 
To give increased evidence to this picture, we have checked numerically the two 
cases next to the M^^^ + 02,1 model, namely M3 7 + 0i,5 and M4J + 02, i- In the first 
case the perturbing operator is even, the scaling function has an expansion both in 
even and odd powers of A and it is not difficult to check that passing from massive 
to massless behaviours amounts to sending A -^ —A. This is common to all p even 
models. In the p odd cases we are in the same situation examined earlier in the 
paper: it is the A — > iA transformation that leads from massive (imaginary A) to 
massless (real A) behaviour. 

5.1 Non-perturbative nature of the flows 

One could ask why we do not check the existence of these massless flows in a 
perturbative framework at large p, like it has been done in unitary minimal models 



perurbed by 01 3 |^, |3^. The reason is that perturbative calculation always fails 
to pick up the IR fixed point in our case, as it turns out from the following analysis. 
The perturbative expansion of the two-point function {ip{x)ip{Q)) around the 
conformal point is given by 

(^(x)V'(O)) = (^(x)V^(0))cFT + A/A(^(a:)^(0)^(2/))cFT (47) 



+ iA2y'rfV'«^(^(a;)V^(0)^(y)V^(t/;))cFT + ... (48) 

Inserting this in the Callan-Symanzik equation provides a perturbative definition of 
the anomalous dimension of the field il){x) and hence the perturbative /5-function [Rl, 



M- 



Consider first the models with p even: ip = 0i 5. For p ^ 00 the parameter 
controlling the perturbation theory is e = 1 — Ai^s = -|^. The first nontrivial 
contribution comes from the order A in the expansion (^8]). The integration is done 
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explicitly in this case leading to the following expression for the /^-function 

P[g)=eg-TxCg^ + 0{g') (49) 

C is the structure constant of the channel 0i,50i,5 — ^ (^01,5 and g the renormalized 
coupling constant. Surprisingly, unlike the case of unitary minimal models per- 
turbed by 01^3, one finds here that C = \e + Oi^e^). We are interested in finding 
a non-trivial fixed point g* in the perturbative region around (7 = 0, i.e. g* ~ e. 
However, formal inserting of the value for C in (l49|) leads to a "nontrivial fixed 
point" 

(?* = ^7r~0(l) (50) 

which is obviously out of the range of perturbation theory. 

Let us now turn to the second case p odd {ip = 02, i)- As it was mentioned above, 
the operator 02_i is Z2-odd and hence its 3-point function vanishes identically. Thus, 
the first nontrivial contibution here comes from the second order, i.e. the 4-point 
function. We analyse the leading term of the channel 02,i02,i —^ -D03,i. This gives 

h = -A2D2(a;2)-2A+^ J d'w\w\^(-''+l')\l - t^l^l^'H J d'y\y - wx\'^''+'^ (51) 

The integration gives a result of order 0(1) and for the structure constant one finds 
explicitely D = \/3 + 0{e). With this result for the two-point function (^Sf ) we shall 
have the following expansion for the /3-function 

P{g)=eg + eAg' + 0{g') (52) 

where A is a numerical coefficient of order 0(1). Again as before a naive compu- 
tation leads to a "fixed point" g* ~ 0(1). We have to conclude that if the models 
considered above actually describe a RG fiow from Mp to Mp_i as indicated by the 
TEA analysis, it should have an essentially nonperturbative nature. 

One remark is in order. The difference in the central charges of the models Mp 
and Mp_i connected by the RG trajectory is 

Ac = — -^ ■—- ~ — ~ 0[e] when p — >• oo (53) 

p{p^ — 1) p 
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At the same time, as we argued above, there is no nontrivial fixed point pertur- 
batively near the origin g = 0. One possible explanation of this fact is as follows. 
The c-theorem of Zamolodchikov is no more valid in our case of nonunitary models. 
But if we fix the "metric" G in the space of the coupling constants to a (positive or 
negative) number by a suitable choice of basis, then 

Ac = -12G r (3{g)dg (54) 

Jo 

Inserting the explicit expression for i3{g) ( ^9]) or (|52| ) one can convince himself that 
it is exactly the value g* ~ 0(1) that ensures Ac ~ 0{e)\ 

6 Conclusions speculations and generalizations 

In this paper we have learned that, allowing the perturbing parameter to take 
complex values, one can greatly enrich the phenomenology of the RG space of 
actions of two-dimensional Quantum Field Theory. This goes in the same direction 
of some recent work of Fendley, Saleur and Al. Zamolodchikov |^, |^ on the Sine- 
Gordon model with imaginary coupling (mass). In the case examined in the present 
paper we are lucky enough to be able to find separately reasonable TBA systems 
for the description of the massive and massless regimes of the models. This allows 
the very detailed analysis done throughout the paper. In other cases, including the 
Sine-Gordon one of [Q, this could not be so easy to imagine at first glance, and one 
has to resort to numerical analytic continuation (Fade approximants) to explore 
the possibilities in the complex A plane. This gives in general less accurate results, 
however in most situations the precision for the IR central charge is still good 
enough to establish some interesting results, as in parafermionic theories perturbed 
by their generating parafermion PB|. 



The possibility to have a massless fiow in a certain direction in the complex 
A plane with a Z2 odd perturbing field is restricted to the cases with conformal 
dimension A = A > |. This guarantees that the condition found in [^ to have 
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square root singularities in the negative r axes does not hold, thus making it possible 
to have a singularity free flow. Notice that the condition A = A > ^ holds neither 
for the 012 nor for the 02i perturbations in the minimal unitary series. We have 
therefore to disappoint our readers about exciting possibilities like critical Ising 
model in imaginary magnetic field, etc... Possible candidates for massless flows 
generated by 02i are the Mp^q models with the condition 

6p > 3g > 4p (55) 

It should be interesting to explore which of them really admit a massless flow for 

Other interesting candidates for this phenomenon can be found outside the range 
of minimal models. Substituting the Ti x A^ structure of TBA with a T„ x A^, one 
could speculate about similar series of massive-massless flows appearing for each 
Wn+i organized non-unitary minimal models. All these series should be related to 
A2n affine Toda QFT, much the like the Mp + ijj are related to Izergin-Korepin. 
Also, many unitary theories perturbed by a relevant Z2 odd field with A > | 
could be considered for this phenomenon. Among these, the most interesting are, 
perhaps, the ^'^^^'^ models, perturbed by their 0^j* operator. The list of possible 
generalizations is very long. 

An important direction of investigation that should be developed is the obser- 
vation that, in analogy with what has been done in 0], undoing the truncation in 
Mp + ip could lead to c-theorem violating flows between different points of some 
c = —2 critical line, interpreted as a kind of "Izergin-Korepin at imaginary cou- 
pling" , exactly as reported in Q for the Sine-Gordon model. The possible relevance 
of these flows for random walks and polymer physics should be considered with some 
attention. 

Another interesting issue is the violation of the Zamolodchikov c-theorem in 
many of these flows. The function c(t) (t a RG "time" ) that decreases monotonically 
in the c-theorem, is not defined in the same manner as the scaling function of 
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TBA. However one can think that the two differ by, say, a "different choice of 
Renormahzation scheme", although this sentence has poor meaning in this context. 
So we expect that the qualitative behaviour of the two are the same, thus the 
conjecture that also the TBA scaling function should decrease monotonically, which 
is confirmed by all the cases of unitary theories analyzed by TBA so far. It is a 
known fact that the c-theorem as it is, holds only for unitary theories, no surprise 
then if half of the models of our Mp + ip series violate it. However people have 
speculated for a long time about some c-theorem, that could hold in non-unitary 
theories if we substitute the scaling function c{r) with the effective one c{r). We 
observe in all our Mp + ip models that such a c-theorem holds. This is true also for 
all the generalizations (like T„ x A^.) where one can explicitly write a massless TBA 
in terms of right and left movers. The cases that violate c-theorem too seem to us 
to be those justified only by analytic continuation, but not admitting a right-left 

mover symmetric TBA, which should also give a signal that 5*+- 7^ S ^. What 

is the deep physical meaning of this fact (parity violation?) should be matter of 
interesting further investigation. 
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